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NON-INVASIVE PIPELINE INSPECTION SYSTEM 

Reference to Related Application 

This application claims priority to the filing date of U.S. Provisional Patent 
Application Serial No. 60/222,170 entitled "Non-Invasive Pipeline Inspection Using 
Radiosounding," filed on August 1, 2000, the disclosure of which is hereby incorporated 
by reference. 

Field of the Invention 

This invention relates generally to inspecting a pipeline for anomalies, and more 
specifically to inspecting a pipeline using a reflected component of an input waveform. 

Background of the Invention 

To maintain substantial fluid flow through a pipeline, internal pipeline 
characteristics need to be monitored so that defects, obstructions, and other anomalies in 
the pipeline can be detected and repaired efficiently. In addition to obstructions affecting 
fluid flow in the pipeline, a pipeline may bend and/or buckle when it experiences a 
change in pressure, such as when the pipeline is laid underwater. Frequently, companies 
must endure substantial monetary costs and schedule delays due to the detection and 
repair of these pipeline anomalies. 

In some conventional systems, an internal, invasive device that crawls the length of 
the pipeline is used to inspect a pipeline for anomalies. This device, often called a "pig", 
poses a serious blockage to the normal fluid flow through the pipeline. A pig may 
additionally require several days for the inspection of a lengthy pipeline. Furthermore, 

1 



the amount of data a pig can record, the life of its battery, and the wear of its components 
from crawling the pipeline all limit the usefulness of the pig. 

Measuring the acoustic signature of a pipeline is another technique used to detect 
pipeline anomalies. This technique sometimes involves hitting the pipeline on its side 

5 with a hard object, such as a hammer, and then measuring the acoustic signature of the 
pipeline. Anomalies often alter the acoustic signature of a pipeline as compared to a 
pipeline with no such anomalies. However, the magnitude of the anomaly that may be 
detected is dependent upon the wavelength of the waveform transmitted along the 
pipeline, and sound waves generally have longer wavelengths than some other 

10 waveforms. Therefore, this technique typically fails to detect smaller-sized anomalies in 
a pipeline. 

Pulse propagation may also be used to detect pipeline anomalies. According to one 
technique, two pulses are transmitted along the pipeline from opposing locations towards 
an intersecting location. The pulses intersect and are each modified by collision with the 

15 oppositely directed pulse. A receiver is positioned at the intersecting location and, after 
receiving the modified pulses, analyzes at least one indicator characteristic of one of the 
modified pulses to determine whether an anomaly exists between the receiver and the 
corresponding transmitter. However, this technique usually requires two separate 
transmitters and a separate receiver, each of which increases the costs associated with 

20 detecting anomalies. Also, pulse propagation analysis may further require inserting the 
receiver into a location in the pipeline not normally open for device placement. 

Another conventional approach is an ultrasonic guided wave inspection technique 
that uses stress waves, such as Lamb waves. Since Lamb waves are typically guided 
along the pipeline, lateral spreading of the energy associated with these waves does not 
25 usually occur and the propagation is essentially one-dimensional. For this reason, Lamb 
waves normally propagate over longer distances than other types of waves, such as bulk 
waves. Unfortunately, at least two modes typically exist at any frequency for Lamb 
waves. Furthermore, the modes are generally dispersive, which means that the shape of 
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the propagating waveform varies with distance along the pipeline. Consequently, 
interpretation of the signals is difficult and can also lead to signal-to-noise problems. 

Accordingly, it is desirable to produce a system that is capable of detecting an 
internal characteristic of a pipeline in a non-invasive fashion. It is also desirable to be 
5 able to inspect a pipeline faster than currently possible, as well as to be able to accurately 
detect smaller-sized anomalies in a pipeline. 

Summary of the Invention 

Briefly, the invention relates to a system and method for inspecting a pipeline. In 

10 one embodiment, the invention provides a system for detecting and characterizing an 
anomaly in a pipeline. In another embodiment, the invention provides a system that can 
also determine the longitudinal path/shape of the pipeline. With a starting point and the 
longitudinal shape of the pipeline, a further embodiment of the invention can also 
determine the location of a pipeline buried underground or even under water. According 

15 to one preferred embodiment, the system includes a processor, an analyzer, and a wave 
launcher. In an alternate embodiment, the analyzer, wave launcher, and processor are 
incorporated into a single unit, thereby eliminating the external connections between the 
devices. The wave launcher communicates with the pipeline, and is adapted to transmit 
an input waveform having a selected input energy along the central longitudinal axis of 

20 the pipeline. Examples of the type of input waveform include, but are not limited to, an 
electromagnetic waveform, a wideband waveform, and an acoustic waveform. Further 
examples of input wideband waveforms include, but are not limited to, a chirp waveform, 
a spread spectrum waveform, a wavelet waveform, and a solitons waveform. The wave 
launcher is further adapted to receive a reflected component of the input waveform 

25 having a characteristic reflected energy. An example of the wave launcher includes an 
antenna adapted to transmit the input waveform along the pipeline. 

The analyzer communicates with the wave launcher, and is adapted to generate the 
input waveform. The analyzer is further adapted to receive the reflected component of 
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the input waveform from the wave launcher. According to a further feature, the analyzer 
includes a signal generator, energy component devices, and a directional coupler. The 
signal generator generates the input waveform that is transmitted along the pipeline. The 
energy component devices extract out the magnitude and phase components of the input 
5 energy and the reflected energy associated with the input waveform and the reflected 
component of the input waveform. An example of the analyzer includes an automated 
vector network analyzer. 

According to one embodiment, the processor communicates with the analyzer, and 
is adapted to compare the input waveform with the reflected component of the input 

10 waveform to generate a mathematical model for the pipeline. According to one feature, 
the mathematical model includes information regarding the longitudinal path/shape of the 
pipeline. According to another embodiment, the processor and the analyzer interact to 
detect and determine the characteristics of an anomaly in the pipeline. Examples of an 
anomaly in the pipeline include, but are not limited to, an obstruction, a flange, rust, and 

15 poorly constructed welds. Specifically, the characteristics of the anomaly include, but are 
not limited to, the location of the anomaly in the pipeline, the type of anomaly, and the 
size of the anomaly. According to a further feature, the system of the invention displays 
the characteristics of the anomaly and/or the shape/location of the pipeline to a user. 

In one operational embodiment, the processor initializes the analyzer. Optionally, 
20 the processor initializes the analyzer by calibrating it. Illustratively, the processor 

calibrates the analyzer by temperature stabilizing the analyzer. In one embodiment, once 
the processor calibrates the analyzer for operation, the user of the inspection system of 
the invention inputs the diameter of the pipeline into the processor. The processor uses 
the diameter to determine the frequency range at which the input waveform can 
25 propagate along the central longitudinal axis of the pipeline. The processor transmits this 
frequency range to the analyzer and the analyzer generates the input waveform having a 
frequency within this range. The analyzer transmits the input waveform to the wave 
launcher, and the wave launcher launches the input waveform along the central 
longitudinal axis of the pipeline. 
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In a further embodiment, the analyzer extracts the input energy associated with the 
input waveform. When the analyzer receives the reflected component of the input 
waveform, the analyzer extracts the reflected energy associated with the reflected 
component. According to another feature, the analyzer then determines a mathematical 
5 representation, or transfer function, relating the input energy and reflected energy. The 
analyzer then transmits this information to the processor. 

In one embodiment, the processor determines the energy reflected from the 
obstruction by generating a mathematical model of the inspection system of the invention 
and the pipeline. According to this feature, the processor determines a model transfer 

10 function relating the input energy of a model input waveform and the model reflected 
energy associated with a model reflected component. In one embodiment, the processor 
determines the characteristics of the obstruction from the model transfer function and the 
transfer function for the input waveform. In some embodiments, the processor 
determines which mathematical model to use from the transfer function relating the input 

15 energy and reflected energy. The processor can use, for example, an ideal physics-based 
model, an average model, and/or a section-by-section model to model the pipeline. 

According to a further feature, the processor displays the characteristics of the 
anomaly with a textual representation on an output device. Alternatively, the processor 
displays the characteristics with a graphical user interface, a three-dimensional solids 
20 rendering plot, or an echo plot. 

In some embodiments of the invention, the pipeline curves along a longitudinal 
central axis. In one aspect, the user of the inspection system still provides the diameter of 
the pipeline. Using the diameter, the processor determines a different range of 
frequencies at which the input waveform can propagate to aid the processor to model the 
25 curved pipeline accurately. Once the processor models the curved pipeline accurately, 
the processor determines the shape of the pipeline, the location of the pipeline, or the 
characteristics of the obstruction. 
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Brief Description of the Drawings 



Advantages of the invention will become better understood by referring to the 
following drawings, which show a system according to an illustrative embodiment of the 
invention and in which: 

5 Figure 1 is a conceptual block diagram depicting the use of a system constructed in 

accord with an illustrative embodiment of the invention for inspecting a pipeline and/or 
for determining the curvature of the pipeline along a longitudinal central axis; 

Figure 2 is a conceptual diagram depicting illustrative waveforms transmitted from 
and received by an exemplary wave launcher used in the system of Figure 1; 

10 Figure 3 depicts an equivalent model of the conceptual diagram of Figure 1 ; 

Figure 4 is an illustrative block diagram showing a lossy physics-based model of 
the conceptual diagram of Figure 1; 

Figure 5 is a flow diagram depicting an illustrative operation of the system of 
Figure 1; 

15 Figure 6 is a more detailed block diagram of the illustrative analyzer of Figure 1; 

Figure 7 depicts a graph describing a probability that a single anomaly will be 
detected using the illustrative system of Figure 1 as the distance between the anomaly and 
the system of Figure 1 increases; and 

Figure 8 is a graph portraying a probability that a single anomaly of varied sizes 
20 (small, medium, large) will be detected using the illustrative system of Figure 1 as the 
distance between the anomaly and the system increases. 

Description of the Illustrative Embodiments 

Figure 1 is a conceptual block diagram depicting an illustrative system 100 for 
inspecting characteristics of a pipeline 140. The illustrative inspection system 100 
25 includes a processor 1 10, an analyzer 120, and a wave launcher 130. In another 
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embodiment, the processor 1 10 is incorporated within the analyzer 120, thereby 
eliminating the external connection between the two devices. In a further embodiment, 
the processor 1 10, the analyzer 120, and the wave launcher 130 are incorporated within a 
single device. As skilled artisans will appreciate, various components of the inspection 

5 system 100 can be implemented in hardware, software or both. The pipeline 140 is 
included in Figure 1 for clarity, but it is not a component of the illustrative inspection 
system 100. Preferably, the inner surface of the pipeline 140 is sufficiently conductive to 
support input waveforms and functions as a waveguide for a suitable axial distance along 
the pipeline 140. As skilled artisans will appreciate, a sufficiently conductive material 

10 may be any of one of a variety of materials, such as, but not limited to, iron, steel, cobalt, 
nickel, and alloys thereof. 

As discussed more fully below, the wave launcher 130 transmits an input 
waveform along the central longitudinal axis 142 of the pipeline 140. In one 
embodiment, the wave launcher 130 is an antenna. The analyzer 120, which is in 

15 communication with the wave launcher 130, generates an input waveform and transmits 
it to the wave launcher 130. According to the illustrative embodiment, the input 
waveform is a wideband waveform, which is a waveform having a bandwidth that forms 
at least about 0.1% of its center frequency. An example of a wideband waveform is a 
waveform which distributes its energy uniformly between 250 MHz and 750 MHz, 

20 having a ratio of bandwidth to center frequency equal to 1 .0 [(750 - 250) / 500 = 1 .0)]. 
In exemplary embodiments, the system 100 employs input waveforms having a center 
frequency of about 800 MHz and a bandwidth of about 400 MHz. Examples of potential 
input waveforms include, but are not limited to, electromagnetic and acoustic waveforms. 

The processor 110, which is in communication with the analyzer 120, processes 
25 and outputs results of the inspection of the pipeline 140. According to one illustrative 
embodiment, the characteristic to be detected is the curvature of the pipeline 140 along 
the longitudinal central axis 142. According to another illustrative embodiment, the 
characteristic to be detected is the diameter of the pipeline 140. In a further illustrative 
embodiment, the characteristic to be detected is the shape of a cross-sectional view of the 
30 pipeline 140, taken for example along view 144. According to the illustrative 
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embodiment of Figure 1, the characteristic of the pipeline 140 to be detected is an 
anomaly in the pipeline 140. In one embodiment, the anomaly is an obstruction 150. In 
other embodiments, the anomaly may be a flange, rust, partly constructed welds, or the 
like. In operation the inspection system 100 detects an illustrative obstruction 150 of the 
5 pipeline 140 that is located a distance 160 away from the wave launcher 130. In 
operation, the analyzer 120 selects the amount of input energy to transmit along the 
pipeline 140. 

In an alternate embodiment, the inspection system 100 determines from one known 
point the location of any other point along the pipeline 140. According to another 
10 embodiment, the inspection system 100 determines the shape (i.e., curvature) of the 
pipeline 140. 

Figure 2 is a conceptual diagram 200 depicting an illustrative input waveform 235 
transmitted from the wave launcher 130, along with an exemplary reflected component 
245. As depicted in Figure 2, the analyzer 120 generates the input waveform 235 

15 corresponding to a selected input energy. The analyzer 120 transmits the input waveform 
235 to the wave launcher 130, and the wave launcher 130 then launches the input 
waveform 235 along the longitudinal central axis 142 of the pipeline 140. After sending 
the input waveform 235 into the pipeline 140, the wave launcher 130 receives a reflected 
component 245 of the input waveform 235. The reflected component 245 includes a 

20 reflected component 245A and a reflected component 245B. The reflected component 
245 A is the component of the input waveform 235 that the obstruction 150 reflects 
towards the wave launcher 130. The reflected component 245B is the component of the 
input waveform 235 that the end wall of the pipeline 140 reflects towards the wave 
launcher 130. The reflected component 245 of the input waveform 235 has a 

25 characteristic reflected energy that depends on the characteristics of the obstruction 150, 
the characteristics of the pipeline 140, the distance 160 between the wave launcher 130 
and the obstruction 150, and other attributes of the illustrative inspection system 100 and 
pipeline 140. These dependencies are further described below. 
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Once the wave launcher 130 receives the reflected component 245 of the input 
waveform 235, the wave launcher 130 transfers it to the analyzer 120, The analyzer 120 
determines the characteristic reflected energy of the reflected component 245 and 
transmits the reflected energy and the input energy to the processor 110. The processor 
5 110 compares the input energy and reflected energy to determine the attributes of the 
obstruction 150. The attributes of the obstruction 150 may be any one of a variety of 
attributes, such as, but not limited to, the size of the obstruction 150, the type of the 
obstruction 150 (e.g., flange, rust, etc.), and the distance 160 to the obstruction 150. The 
processor 110 then reports its results on an output device connected to the processor 1 10. 

10 According to a further feature, the illustrative processor 1 10 begins by calibrating 

the analyzer 120 for measurement. In one embodiment, the processor 110 calibrates the 
analyzer 120 by temperature stabilizing the analyzer 120. Temperature stabilizing 
includes an operator of the illustrative system 100 positioning the analyzer 120 in a 
temperature cycling chamber. In one embodiment, the temperature cycling chamber is an 

15 enclosed, insulated area that introduces devices such as an analyzer 120 to a range of 
temperatures. The processor 1 10 is positioned outside of the temperature cycling 
chamber. The processor 110 loads from its processor memory (e.g., ROM, RAM) a test 
program at which the analyzer 120 can perform several functions and operations 
specified in the test program. For example, the processor 110 may request the analyzer 

20 120 to perform the operations corresponding to the future operations that the analyzer 
120 will carry out. Alternatively, the processor 110 may request the analyzer 120 to 
perform a diagnostic test on the components of the analyzer 120. 

The processor 1 1 0 begins this test program and subsequently introduces the 
analyzer 120 to a range of temperatures while the analyzer 120 is in operation. Once the 

25 analyzer 120 is subjected to the entire range of temperatures, it becomes temperature 
stabilized and it transmits the results from the test program to the processor 110. The 
processor 110 receives and stores the results of the analyzer 120 running this test 
program. When the illustrative inspection system 100 is later positioned at the pipeline 
140, the processor 110 measures the ambient temperature at the pipeline 140. The 

30 processor 110 consequently retrieves the stored results of the analyzer 120 in the 
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temperature cycling chamber running the test program at the ambient temperature. The 
processor 110 initializes the analyzer 120 by using the stored results at the ambient 
temperature. 

In another embodiment, the processor 110 calibrates the analyzer 120 every time 
5 the analyzer 120 is powered up. As described above, the processor 110 measures the 
ambient temperature of the pipeline 140 and executes the test program on the analyzer 
120. The analyzer 120 executes the test program at the current temperature and then 
transfers the results to the processor 110. The processor 110 compares these results with 
expected results at the ambient temperature to obtain a temperature error associated with 
10 the analyzer 120. The processor 110 calibrates the analyzer 120 in this fashion every 

time the temperature at the location at which the illustrative inspection system 100 is used 
varies from the previous temperature at the previous location. Moreover, the processor 
110 calibrates the analyzer 120 in this fashion whenever the analyzer 120 is powered 
down and then powered up. In a further embodiment, the processor 110 alerts the 
15 operator of the illustrative system 1 00 when the temperature error is above a 
predetermined temperature error threshold. 

According to another embodiment, the processor 110 calibrates the analyzer 120 
by temperature stabilizing the analyzer 120 in a thermostatically-controlled chamber. In 
one embodiment, the thermostatically-controlled chamber is a temperature cycling 

20 chamber, as described above, operating at a continuous, constant temperature. By way of 
example, the thermostatically-controlled chamber operates at 25° Celsius. The operator 
of the illustrative system 100 positions the analyzer 120 in the thermostatically-controlled 
chamber and the inspection system 100 begins normal execution. In a further 
embodiment, the processor 110 compares the output of the analyzer 120 at the constant 

25 temperature with expected results at the same constant temperature to obtain a 

temperature error associated with the analyzer 120. In still a further embodiment, the 
processor 110 displays a warning to the operator of the illustrative system 100 when the 
temperature error is above a predetermined temperature error threshold. Alternatively, 
the processor 110 initializes the analyzer 120 with one of the calibration techniques 
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described above or below when the temperature error is above the predetermined 
threshold. 



Once calibration is complete, the processor 110 instructs the analyzer 120 to 
generate an input waveform 235 that will be transmitted along the pipeline 140. The 

5 analyzer 120 may generate an input waveform 235 using a signal generator. 

Alternatively, the analyzer 120 applies a force to the pipeline 140 to generate a sound 
wave as the input waveform 235. The processor 110 indirectly selects the input energy of 
the input waveform 235 by selecting the frequency of the input waveform 235. Before 
transmitting the input waveform 235 to the wave launcher 130, the analyzer 120 

10 determines the input energy associated with the input waveform 235. 

As discussed in more detail below with respect to Figure 6, after the analyzer 120 
determines the input energy for to the input waveform 235, the analyzer 120 transmits the 
input waveform 235 to the wave launcher 130. The wave launcher 130 in turn launches 
the input waveform 235 along the central axis 142 of the pipeline 140. Then the wave 
15 launcher 130 receives the reflected component 245 of the input waveform 235 and sends 
it to the analyzer 120. 

Once the analyzer 120 receives the reflected component 245, the analyzer 120 
determines a transfer function relating the input energy corresponding to the input 
waveform 235 with the reflected energy corresponding to the reflected component 245 of 
20 the input waveform 235. The analyzer 120 determines a transfer function for each 
reflected component 245 (e.g., reflected component 245A and 245B) of the input 
waveform 235. The transfer function of energy is denoted by the following equation: 

transfer function = ^fleeted ^ 

^ input 

Once the analyzer 120 determines a transfer function for the input energy and the 
25 reflected energy corresponding to the reflected components 245A and 245B, it transmits 
these transfer functions to the processor 110. 
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Figure 3 depicts an equivalent model 300 of the illustrative inspection system 100 
and the pipeline 140 of Figure 1. The processor 110 determines the energy reflected 
from the obstruction 150 by generating a mathematically modeled pipeline that is 
representative of the pipeline 140. The analyzer 120 simulates the input waveform 235 

5 that is transmitted along the pipeline 140 as a model input waveform 305. The model 
input waveform 305 is shown at the lower left corner of Figure 3. The analyzer 120 
transmits the model input waveform 305 to the wave launcher 130 in preparation for the 
launching of the model input waveform 305 along the central axis of the model pipeline. 
As a result of imperfections in test port cables and other calibration effects, a calibration 

10 component 3 1 0 of the model input waveform 305 is immediately reflected back to the 
analyzer 120. This calibration component 310 and the energy associated with this 
calibration component 310 is represented in Figure 3 as Hcaiibration(f)- 

A first remainder 320 and a second remainder 330 of the model input waveform 
305 are transmitted through the wave launcher 130 and travel the distance 160 to the 
15 model obstruction 335, or model target. The remainders 320 and 330 are represented in 
Figure 3 as Hiauncherd) and Hpj(f f dj), respectively. The wave launcher 130 has intrinsic 
losses associated with it, and so when the model input waveform 305 is transmitted 
through the wave launcher 130 into the model pipeline, a reflected wave launcher portion 
370 of the model input waveform 305 gets reflected toward the analyzer 120. 

20 At the distance 160, the model obstruction 335 causes a first model reflected 

component 333 of the model input waveform 305 to be reflected toward the wave 
launcher 130. The first model reflected component 333 represents the reflected 
component 245 A shown in Figure 2. A third remainder 340 of the model input waveform 
305 continues along the model pipeline until it reaches the end of the model pipeline. A 

25 second model reflected component 350 is then reflected toward the wave launcher 130 
when it reaches the end wall of the model pipeline, and this second model reflected 
component 350 represents the reflected component 245B. The second model reflected 
component 350 and the energy corresponding to this reflected component 350 is 
represented in Figure 3 by Hp2(f, d^). The sum of the model reflected components 33 3 , 

30 350 are combined at a first summation block 355 and the resulting sum 360 is transmitted 
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to the analyzer 120. The resulting sum 360 is transmitted through the wave launcher 130, 
which is shown in Figure 3 as H P j(f t dj). Additionally, the model obstruction 335 reflects 
a portion of the second model reflected component 350 (that was reflected by the end 
wall of the model pipeline) back toward the end wall, creating a third model reflected 
5 component 353. 

The wave launcher 130 transmits the resulting sum 360 and the reflected wave 
launcher portion 370 toward the analyzer 120. The resulting sum 360 and the reflected 
wave launcher portion 370 are combined with the calibration component 310 and any 
analyzer 120 noise sources 315 at a second summation block 375. A total model 

10 reflected component 380 is then transmitted to the analyzer 120. The total model 

reflected component 380 therefore includes a reflected component corresponding to the 
wave launcher 130 (e.g., reflected wave launcher portion 370), the model obstruction 335 
(e.g., first model reflected component 333), the end wall of the model pipeline (e.g., 
second model reflected component 350), the calibration effects (e.g., calibration 

15 component 310), and any noise associated with the analyzer 120 (e.g., analyzer 120 noise 
sources 315). 

As described in more detail below with respect to Figure 6, the analyzer 120 
receives the total model reflected component 380 and calculates a model transfer function 
relating the model input energy with the model reflected energy corresponding to the 

20 total model reflected component 380. The analyzer 120 then transfers this model transfer 
function to the processor 110. The processor 110 compares the transfer function 
associated with the reflected energy of the reflected component 245 to the model transfer 
function corresponding to the total model reflected component 380. From this 
comparison, the processor 110 determines the location 160 and size of the obstruction 

25 150 and reports these results on an output device. 

In one embodiment, the processor 110 includes the calibration component 310 of 
the analyzer 120, the response of the wave launcher 130, and the response of the pipeline 
140 stored in its local memory (e.g., RAM, ROM). The analyzer 120 noise may be 
negligible if the pipeline 140 reflects most of the input waveform 235. In this situation, 
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the processor 110 can detect a minute obstruction 150 at virtually unlimited range. In 
another embodiment, the processor 1 10 accounts for the analyzer 120 noise when the 
analyzer 120 receives the reflected component 245. 

Illustratively, to process the calculations and modeling as described above, the 
5 processor 110 has digital signal processing capabilities that are used in a collection of 
DSP algorithms (discussed in further detail below). In one embodiment, the processor 
110 uses an ideal lossless physics-based model as the hypothetical model to represent a 
pipeline 140 with no contaminants, defects, anomalies, or other losses. The model 
pipeline has uniform quality of construction material, an identical cross-section along the 
10 entire length of the model pipeline, and a perfectly conductive inner surface. Because the 
ideal physics-based model includes a model pipeline that is an ideal pipeline 140 with no 
losses, the processor 110 determines the response of the model pipeline and subsequently 
determines the type of the obstruction 150 and the location 160 of the obstruction 150 
within the pipeline 140 by comparing the actual reflected energy of the pipeline 140 with 
15 the modeled reflected energy of the ideal pipeline 140. In another embodiment, the 
processor 110 uses an ideal lossy physics-based model. The processor 110 models the 
pipeline 140 using an ideal lossy physics-based model when the processor 110 assumes a 
conductive inner surface that experiences greater losses relative to the conductivity of the 
inner surface of the model pipeline. 

20 Figure 4 is an illustrative block diagram showing a lossy physics-based system 

model 400 of the inspection system 100 incorporating partial a priori knowledge. As 
previously described, the analyzer 120 generates a series N of input waveforms 235 and 
applies these input waveforms 235 to the wave launcher 130. The amplitude 
x(f n ),n = 0,1 V ..,N - 1 , of the input waveforms 235 is a function of the excitation 

25 frequency of the input waveform 235. The model 400 also includes the reflection 

response 410 of the wave launcher 130 and other near-field effects (i.e., the effects on the 
electric and magnetic fields of the reflected component 245 when the reflected 
component 245 is within the range of the wave launcher 130), denoted below by H B (fn). 
The processor 110 represents partial knowledge of an as-built inspection system 100 in 

30 the reflection response 410 of the wave launcher 130. When reflected toward the wave 



-14- 



launcher 130, the input waveform 235 further experiences a scaling coefficient 420 for 
near-field effects, represented below by Kb. The scaling coefficient 420 adjusts the 
magnitude and phase of the reflected component 245. 

The processor 110 models the pipeline 140 as a lossy physics-based model 425 , 
5 shown as H T (f n ; d, a, a). The lossy physics-based model 425 of the pipeline 140 depends 
on several parameters of the pipeline 140, such as, but not limited to, the round-trip 
distance d between the obstruction 150 and the wave launcher 130, the radius a of the 
pipeline 140, the effective conductance a of the pipeline's 140 inner surface, the scaling 
coefficient 430 Kr for the obstruction 150, and the background noise 435 7](f^ of the 
10 analyzer 120. 

For each of these input waveforms 235, the wave launcher 130 receives the 
reflected component 245 of each input waveform 235. The amplitude y(fy) of the 
reflected component 245 is also a function of the excitation frequency of the input 
waveform 235. The analyzer 120 calculates an estimate of the inspection system 100 
15 transfer function, which, as described above, is given as: 

H{f n )-^y (1) 
Afn) 

The processor 110 then operates on this transfer function to locate and identify any 
obstructions 150 within the pipeline 140. Under the assumption that the background 
noise rjtf^ is a zero-mean, independent complex Gaussian process, the processor 110 
20 employs the minimum mean-squared error estimate, given as: 

J* =Y\H{f^k B H B (fyk T H T {f n] la^ . (2) 

«=0 



Note that, following standard convention, the carat ( A ) calls our attention to an estimated 
value of a parameter (as opposed to its "true" value). 
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To begin signal processing, the processor 110 assumes a range of distances over 
which to search for obstructions 150 within the pipeline 140. This range is denoted as <4 
1 = 0, 1 , . . . , L- 1 , where L is the total number of steps within the range d\ of distances at 
which to search for obstructions 1 50. In one embodiment, the range d\ may cover a few 
kilometers in steps of 0.1 meters. For each value of <4 the pipeline 140 transmission is 
calculated as: 



H T {f n ;d„a,cr) = e- a » d <e- jM , 



(3) 



where 
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15 



In Equations (4) and (5), the new symbols are identified as: 



S 0 Permittivity of free space, 8.854 x 10 12 



ju 0 Permeability of free space, 4tt x 1 0 



N-m 2 
7 Wb 



A • m 

First root of the first derivative of the Bessel function of first kind. 



Given the three vectors, H{f n ), H B (/„ ), H T (f n ;d,,a,&), the processor 1 1 0 
20 calculates the inter- and intra-signal correlation functions as: 
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K B „ T =f,K(fMf n -,d,,a,a), 



(6) 
(7) 
(8) 
(9) 
(10) 



and the measurement correlation functions as: 



10 



P HT =f j H;(f n ;d l ,a,a)Hif ii ). 



(11) 
(12) 



The signal correlation functions are used to form the correlation matrix 



Rrj rr R 



R R 



15 



while the measurement correlation functions are incorporated into the vector 



p(<0= 



(13) 



(14) 



For a particular selection of distance d l , the minimum mean-squared error is given as: 



20 
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(15) 



The associated values of the optimum scaling constants are given as 



K T (d,\ 



R _1 p • (16) 



Equations (15) and (16) are computed for all d n l~ 0,1,..., L-l. Once completed, the 
global minimum attained by J m i n is identified together with the distance d l at which it 
occurs, and the attendant value of Kt. 

The magnitude of the estimate of Kris related to the cross-sectional area of the 
1 0 obstruction 1 5 0 as : 

|lU«8.3x7 ,2 +0.5x7\ (17) 



where T is the fractional cross-sectional area of the obstruction 150. This expression is 
15 inverted to find the size of the target 150. 

To find the type of the obstruction 1 50, the magnitude of J min at its global 
minimum helps define the depth (i.e., distance 160 along the length of the pipeline 140) 
of the obstruction 150. For example, since the processor 1 10 bases the lossy physics- 
based model 425 on the presumption of an obstruction 150 having "zero thickness," an 

20 obstruction 1 50 of substantial length provides a relatively high value at the local 

minimum. In one embodiment, the processor 110 displays the magnitude of J min at its 
global minimum to the user of the system 100 so that the user can determine the type of 
the obstruction 150. In another embodiment, the processor 110 has a table stored in local 
memory associating a range of depths with a type of obstruction 150 and determines the 

25 type of obstruction 1 50 from the depth and the stored table. 

Since J m in is calculated as a function of distance (see Equation 15), the processor 
1 1 0 determines the location 1 60 of the obstruction 1 50 from the distance di. The distance 
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di at which the global minimum of i/min occurs is the maximum likelihood estimate of the 
target 150 location 160 from the wave launcher 130. 

In an alternate embodiment, the processor 110 employs the method of maximum 
likelihood, which requires full knowledge of the output probability density functions, to 
5 locate and identify any and all obstructions 150 within the pipeline 140. The description 
provided above is strictly illustrative and does not make pretense to describe the myriad 
of improvements possible or the field of alternatives available. 

As skilled artisans will appreciate, the processor 110 may be any one of a variety 
of devices, such as, but not limited to, a laptop computer with digital signal capabilities, a 
10 personal digital assistant with digital signal capabilities, a mobile telephone with digital 
signal capabilities, and a beeper with digital signal capabilities. Generally, the processor 
110 can be any device that has computer memory (e.g., RAM, ROM) and digital signal 
processing capabilities so that the DSP algorithms can be stored and executed on the 
processor 110. 

15 In another embodiment, the processor 110 uses an average model of the entire 

pipeline 140. The processor 110 averages losses associated with construction, internal 
characteristics, differences within the cross section, and other losses apparent throughout 
the pipeline 140 to obtain an average model pipeline. In yet another embodiment, the 
processor 110 utilizes a section by section model of the pipeline 140, in which the 

20 processor 110 segments the pipeline 140 into sections and computes a representation for 
each of the segmented sections. The processor 110 builds a model of a portion of the 
pipeline 140 being tested by analyzing and then joining each section of the relevant 
portion of the pipeline 140. 

According to one embodiment, the operator of the inspection system 100 selects 
25 the appropriate model (e.g., ideal physics-based system model 400, average model, 

section-by-section model) that the processor 1 10 uses to model the pipeline 140 from a 
menu displayed on the output device, as discussed more fully below. According to 
another embodiment, the processor 110 determines which model to apply depending on 
the characteristics of the pipeline 140. As described in more detail below with respect to 
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Figure 6, after the processor 110 receives the transfer function relating the input energy 
corresponding to the input waveform 235 with the reflected energy corresponding to the 
reflected component 245 , the processor 110 uses this data to determine which model to 
use in determining the characteristics of the obstruction 150. 

5 Figure 5 is a flow diagram 500 depicting an illustrative operation of the inspection 

system 100 of Figure 1. First, the processor 110 initializes (Step 510) the analyzer 120. 
In one embodiment and as described above, the processor 110 calibrates the analyzer 120 
by temperature stabilizing the analyzer 120. The processor 110 may also initialize the 
analyzer 120 by additionally or exclusively performing a diagnostic check on the 

10 components of the analyzer 120, which are described further in Figure 6. In another 
embodiment, the processor 1 10 initializes the analyzer 120 by supplying power to the 
analyzer 120. Initialization may also include a combination of the techniques described 
above. 

At step 520, the wave launcher 130 transmits the input waveform 235 along the 
15 central longitudinal axis 142 of the pipeline 140. As described above, the analyzer 120 
generates the input waveform 235 and transmits it to the wave launcher 130. More 
specifically, in one embodiment the input waveform 235 is an electromagnetic waveform 
having a selected frequency and energy. The range of frequencies at which the input 
waveform 235 is generated is discussed more fully below with respect to Figure 6. As 
20 skilled artisans will appreciate and as described more fully below, the input waveform 
235 may be any one of a variety of wideband waveforms, such as, but not limited to, a 
chirp waveform, a spread spectrum waveform, a wavelet waveform, and a solitons 
waveform. In another embodiment, the input waveform 235 is not an electromagnetic 
waveform but rather is an acoustic waveform. 

25 After the wave launcher 130 transmits the input waveform 235 along the central 

longitudinal axis 142 of the pipeline 140, the wave launcher 130 receives the reflected 
component 245 of the input waveform 235 and transmits it to the analyzer 120. As 
described above with respect to Figure 2, the analyzer 120 measures (Step 530) the 
characteristic reflected energy of the reflected component 245. The analyzer 120 then 
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determines the transfer function relating the input energy corresponding to the input 
waveform 235 with the reflected energy corresponding to the reflected component 245 
(e.g., reflected components 245A and 245B) of the input waveform 235. The analyzer 
120 then transmits the transfer function to the processor 1 10. 

5 As discussed with respect to Figure 3, the processor 110 compares (Step 540) the 

transfer function associated with the reflected energy of the reflected component 245 with 
the model transfer function corresponding to the total model reflected component 380. 
From this comparison, the processor 1 10 determines (Step 550) the location 160 and size 
of the obstruction 150. Although the flow diagram 500 illustrates the operation of the 

10 inspection system 100 for one obstruction 150, the invention extends to a pipeline 140 
containing a plurality of obstructions 150. In other embodiments, step 550 may also 
include using the above discussed mathematical process to determine the axial shape of 
the pipeline 140 (i.e. the curvature of the pipeline 140 along the central longitudinal axis 
142. Skilled artisans will appreciate that the pipeline 140 need not have a circular cross- 

15 section 144 and that the location of the central axis 142 along which the input waveform 
235 propagates may be adjusted to accommodate pipelines 140 having non-circular 
cross-sections 144. As previously mentioned, in some embodiments, a user provides the 
measurement system 100 of the invention with cross-sectional information of the pipeline 
140. However, according to other embodiments, the system 100 automatically 

20 determines the cross-section 144 of the pipeline 140. 

At step 560, the processor 110 displays the results on an output device to the user 
of the inspection system 100. The reported results may be any of one of a variety of 
statistics, such as, but not limited to, the type of the obstruction 150, the size of the 
obstruction 150, and the location 160 of the obstruction 150. Examples of output devices 
25 are any one of a variety of devices such as, but not limited to, a computer monitor, a LCD 
screen, one or several light sources having a predefined meaning associated with the 
obstruction 150 (e.g., a blue light denoting that the obstruction 150 is a flange, a red light 
indicating that the obstruction 150 is rust, etc.), a cellular phone screen, a personal digital 
assistant screen, and an output device that generates predefined tones (e.g., a 40 Hz tone 
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meaning the obstruction 150 is a flange, a 60 Hz tone meaning the obstruction 150 is rust, 
etc.). 



In one embodiment, the processor 110 displays the results corresponding to the 
obstruction 150 with a graphical user interface (GUI). The output device associated with 
5 the processor 1 10 displays the GUI, and the GUI may represent the obstruction 1 50 with 
images, buttons, scales, etc. Alternatively, the processor 110 displays the results with an 
echo plot, which is a plot that displays points to trace the location 160 and size of the 
obstruction 150 in the pipeline 140. In yet another embodiment, the processor 110 
displays the results with a textual description. For example, the processor 110 reports 
10 that the obstruction 150 is a "3 cm buckle found at 10 km". The processor 110 may also 
report the results with a 3-dimensional solids rendering plot. Although several 
techniques to output the results are described above, skilled artisans will realize that other 
output methods may be used in place of or in combination with the above techniques. 

Figure 6 is a more detailed block diagram 600 of the illustrative analyzer 120 of 
15 Figure 1 . In one embodiment, the analyzer 120 is an automated vector network analyzer. 
The analyzer 120 includes a signal generator 610, energy component devices 620A and 
620B, and a directional coupler 630. The directional coupler 630 transmits an input 
energy 615 associated with the input waveform 235 to the energy component device 
620 A. The directional coupler 630 transmits a reflected energy 625 associated with the 
20 reflected component 245 to the energy component device 620B. The directional coupler 
630 transmits the input waveforms 235 to the wave launcher 130 along a first 
communication channel 635. The wave launcher 130 transmits the reflected component 
245 of the input waveform 235 to the analyzer 120, and more specifically to the 
directional coupler 630, along a second communication channel 640. 

25 The signal generator 610 generates the input waveform 235 that is transmitted 

along the pipeline 140. The signal generator 610 generates input waveforms 235 having 
a frequency within a certain range of frequencies, determined by the characteristics of the 
signal generator 610 and by the characteristics of the pipeline 140. As described above, 
the pipeline 140 acts as a waveguide for the input waveform 235, and input 
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electromagnetic waveforms 235 propagate along a waveguide with different field 
configurations (e.g., electric field and magnetic field) and different velocities. This is 
referred to as the mode of the wave, and different modes of a wave can propagate along a 
waveguide simultaneously. 

5 The pipeline 140 has a cutoff frequency below which no input waveform 235 will 

propagate. This cutoff frequency is the minimum frequency needed to propagate the first 
mode of the input waveform 235 along the pipeline 140. The first mode of an 
electromagnetic waveform, which propagates along the pipeline 140 alone, is called the 
dominant mode of the waveguide. The minimum frequency at which the dominant mode 

10 exists, which is the cutoff frequency, depends on the cross-section 144 of the opening of 
the pipeline 140. The maximum frequency at which the dominant mode exists depends 
on the characteristics of the pipeline 140. 

In one embodiment, the pipeline 140 is a circular cylindrical pipeline 140, and the 
range of frequencies at which the dominant mode propagates is given by the following 
15 relationship: 

1 < f d < 2 



a a 

wherein: 

f d is the frequency at which the dominant mode propagates along the pipeline 140; 

c is the speed of light ( 2.998 x 1 0 s meters / second); 

20 Ki and K2 are constants associated with the characteristics of the pipeline 140; and 

a is the radius of the circular cross-section 144 of the pipeline 140. 

For a circular cylindrical pipeline 140, the dominant mode is referred to as the TEn 
mode. TE waves are waves in which the longitudinal components of the electric field at 
the walls of the waveguide are zero and the longitudinal magnetic field is non-zero. In 



-23 - 



one embodiment, the signal generator 610 transmits the dominant mode of the input 
waveform 235. The illustrative signal generator 610 generates an input waveform 235 
for the entire range of frequencies at which the dominant mode exists. When the input 
waveform 235 is at a frequency associated with the dominant mode, the analyzer 120 
5 generates a unique transfer function relating the input energy 615 and the reflected 

energy 625. The transfer function is unique because the dominant mode is the only mode 
of the input waveform 235 that propagates along the pipeline 140. 

According to one illustrative embodiment, the user of the inspection system 100 
enters the diameter information for the pipeline 140 into the processor 110. According to 

10 another embodiment, the user enters the shape and dimensions of the cross-section 144 of 
the pipeline 140 into the processor 110. The processor 110 uses the entered information 
to determine the frequency range at which the dominant mode of the input waveform 235 
propagates. The processor 110 then notifies the analyzer 120 to generate input 
waveforms 235 each having a frequency within the range of frequencies of the dominant 

15 mode. Alternatively, the user of the inspection system 100 enters the brand name of the 
pipeline 140 and the processor 110 uses this data to retrieve from its local memory the 
cross-sectional information of the pipeline 140. Generally, the user of the inspection 
system 100 can input any parameter of the pipeline 140 into the processor 1 10 as long as 
the processor 1 1 0 can determine the frequency range of the dominant mode of the 

20 pipeline 140. 

In one embodiment and as briefly described above with respect to Figure 5, the 
analyzer 120 generates a chirp waveform as the input waveform 235. A chirp waveform 
is a quasi-sinusoidal waveform that has the property that its instantaneous frequency is a 
linear function of time. The analyzer 120 generates discrete chirp waveforms and 
25 increments the frequency of the input waveform 235 by a step-size through a range of 
sinusoidal frequencies. By way of example, the analyzer 120 generates discrete chirp 
waveforms and increments the frequency by a step-size of 1 Hz through 3 Hz (i.e., the 
analyzer 120 generates discrete chirp waveforms having a frequency of 600 MHz, 601 
MHz, and 602 MHz). 
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In another embodiment and as briefly described above with respect to Figure 5, the 
analyzer 120 generates a prototype waveform and derives a wavelet waveform as the 
input waveform 235. The analyzer 120 derives a wavelet waveform by stretching or 
delaying the prototype waveform. The analyzer 120 has a high degree of control over the 
5 joint time and frequency distribution of the input energy 61 5 in the wavelet waveform. 
For example, a wavelet waveform can be derived such that all frequency components 
arrive at substantially the same time and substantially in phase. 

In another embodiment, the analyzer 120 generates a spread spectrum waveform as 
the input waveform 235. The spread spectrum waveform reduces interference by 

10 spreading the input waveform 235 in bandwidth prior to transmission along the pipeline 
140. Upon receiving the reflected component 245 of the input waveform 235, the 
analyzer 120 despreads, or decreases, the bandwidth of the reflected component 245 by 
the same amount of bandwidth as the increase. This technique in turn decreases the 
effect of the interference that occurs during the transmission and reception of the input 

15 waveform 235 and the reflected component 245. 

When the wave launcher 130 launches many input waveforms 235 of different 
frequencies within the range of frequencies of the dominant mode, each input waveform 
235 travels along the central axis 142 of the pipeline 140 at different velocities due to the 
different frequencies. This is referred to as "dispersion" of the input waveform 235. 

20 When the pipeline 140 is a relatively straight pipeline 140, the operation of the inspection 
system 100 is not affected by the different velocities of the input waveforms 235 because 
each input waveform 235 has a separate component 245 of the input waveform 235 
reflected toward the wave launcher 130 at different times corresponding to the different 
velocities. Therefore, the inspection system 100 detects the obstruction 150 when the 

25 input waveform 235 disperses in a straight pipeline 140. 

In another embodiment, the pipeline 140 is a pipeline 140 that has curves and 
bends. As previously described above, the user of the inspection system 100 may provide 
information such as cross-sectional and axial curvature information to the processor 110. 
The processor 110 uses this information to calculate the range of frequencies 
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corresponding to the dominant mode as well as the range of frequencies corresponding to 
higher order modes of the input waveform 235 and to generate a mathematical model of 
the pipeline 140. Alternatively, the system 100 determines the cross-sectional and axial 
curvature properties of the pipeline 140. Either way, in one embodiment the signal 
5 generator 610 generates input waveforms 235 within a range of frequencies that 

correspond to more than one mode of the input waveform 235 (i.e., the dominant mode 
and higher order modes). The wave launcher 130 then launches these input waveforms 
235 along the central axis 142 of the curved pipeline 140. The analyzer 120 receives an 
independent reflected energy 625 along the second communication channel 640 for each 
10 input waveform 235 that was introduced. 

In one embodiment and as described above, the processor 110 compensates for 
dispersion in its formulation of the model pipeline and therefore forces time-alignment of 
all the frequencies of the input waveforms 235 that travel at different velocities. The 
pipeline 140 incorporates dispersion into its DSP algorithms to model the pipeline 140 
15 because the dominant mode dispersion of an input waveform 235 is substantially 

identical in both a straight and curved pipeline 140. For example, the lossy physics-based 
model 425 described above compensates for dispersion. More specifically, the lossy 
physics-based model 425 described above incorporates dispersion in its formulation of 
the model pipeline with the term under the second radical in Equation (5). 

20 In another embodiment, the processor 1 1 0 uses the transfer function of each input 

waveform 235 to determine which model (ideal physics-based system model 400, 
average model, section-by-section model) of the pipeline 140 to use. Therefore, the 
analyzer 120 helps the processor 110 accurately model the pipeline 140 when the 
analyzer 120 generates higher order mode input waveforms 235 for a curved pipeline 

25 140. 

The processor 110 models the curves in a pipeline 140 more realistically as the 
number of modes that are propagating increases because of dispersion, which was 
described above. As the frequency of the input waveforms 235 increases, and therefore 
higher order modes propagate, the input waveforms 235 propagate around curves with 
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greater differences in velocities relative to the difference in velocities along a relatively 
straight portion of the pipeline 140. The processor 1 10 models the curves more 
accurately due to these velocity differences. Therefore, the inspection system 100 detects 
the obstruction 150 when the pipeline 140 is a curved pipeline 140. 

5 In another embodiment and as briefly described above with respect to Figure 5, the 

analyzer 120 generates a soliton waveform as the input waveform 235. A soliton 
waveform is a class of waveforms designed to pass through a non-linear dispersive media 
without losing its shape and properties. The processor 110 uses soliton waveforms as the 
input waveform 235 to characterize the curvature of the pipeline 140. In one illustrative 

10 approach, the processor 1 1 0 determines the curvature of the pipeline 140 by refining the 
shape of the soliton waveform in real-time until the analyzer 120 receives an unchanged 
reflected component 245. Alternatively, the processor 110 refines the spectral content of 
the soliton waveform in real-time until the analyzer 120 receives an unchanged reflected 
component 245. In another embodiment, the processor 110 refines the power level of the 

15 soliton waveform in real-time until the analyzer 120 receives an unchanged reflected 
component 245. 

In another embodiment, the pipeline 140 is a hollow rectangular pipeline 140, and 
the dominant mode of the input waveform 235 propagates over the range of frequencies 
given by the following relationship: 

C r C 

20 < J d 

2a 2b 

wherein: 

f d is the frequency at which the dominant mode will propagate along the rectangular 
pipeline 140; 

c is the speed of light (2.998 x 10 8 meters / second); 
25 a is the height of the pipeline 140; and 
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b is the width of the pipeline 140, assuming the width is less than the height of the 
pipeline 140. 

According to this embodiment, the user of the inspection system 100 provides the 
processor 1 1 0 with the height and width of the pipeline 140. With these parameters, the 
5 processor 110 determines the range of frequencies at which the dominant mode and 
higher order modes of the input waveform 235 propagate along the hollow rectangular 
pipeline 140. 

The energy component devices 620A, 620B extract out the magnitude and phase 
components of the input energy 615 and the reflected energy 625 associated with the 

10 input waveform 235 and the reflected component 245, respectively. The processor 110 
requires the magnitude and phase of the input energy 615 and the reflected energy 625 to 
determine the attributes of the obstruction 150. The energy component devices 620A and 
620B do not affect the input waveform 235, the reflected component 245, the input 
energy 615, or the reflected energy 625 when extracting out the magnitude and phase of 

1 5 the input energy 6 1 5 and the reflected energy 625 . 

The directional coupler 630 transmits and receives energy between the signal 
generator 610, the energy component device 620B, and the wave launcher 130 without 
any physical connection between the devices. In one embodiment, the directional coupler 
630 uses the electric fields generated by the circuits of these components to transmit and 
20 receive energy. 

According to a further embodiment, the inspection system 100 detects the axial 
curvature of the pipeline 140 with or without an obstruction 150. As described above, the 
wave launcher 130 launches input waveforms 235 corresponding to the dominant mode 
and higher order modes of the input waveforms 235 along the central axis 142 of the 
25 pipeline 140. 

The axial curvature of the pipeline 140 may be useful to the user of the inspection 
system 100 for a variety of reasons. By way of example, it can be useful to determine a 
change in the degree of curvature over a period of time and to locate the end wall of the 
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pipeline 140 when the end wall is not located at the expected location, and the like. The 
change in the degree of curvature over a period of time also shows, for instance, a portion 
of the pipeline 140 experiencing a greater amount of force applied to it relative to less 
curved portions. The user can use the curvature information to adjust characteristics of 

5 the pipeline 140 such as re-position the pipeline 140 in a modified location, pad the 

curved portion of the pipeline 140 to adjust (i.e., decrease) the amount of force applied to 
it, apply a similar force to the uncurved portions of the pipeline 140 to decrease the rate 
of curvature change along the pipeline 140, or the like. In yet another embodiment, 
knowing an initial location point along the pipeline 140 and curvature information 

10 determined by the inspection system 100, a user of the inspection system 100 can map the 
location of a length of pipeline 140, even if the pipeline 140 is underground or submerged 
underwater. 

The location of any point along the pipeline 140, such as the point corresponding to 
the end wall of the pipeline 140, may be useful to the user of the inspection system 100 

15 for a variety of reasons. For example, a user of the inspection system 1 00 may know the 
location of a point along the pipeline 140 but may not know the location of the end wall 
of the pipeline 140 if the pipeline 140 is laid underground or underwater. Similarly, 
although a user of the inspection system 100 may know the starting point of an old 
pipeline 140 buried in the foundation of a building, a user may not know the path the 

20 pipeline 140 takes throughout the foundation. One skilled in the art will appreciate that 
knowing the location of an entire segment of pipeline 140 may, for example, aid in repair 
of an anomaly 150. Such information may also be helpful with regard to installing 
additional pipeline 140 segments. 

Figure 7 shows a graph 700 describing the probability that the inspection system 
25 100 detects the obstruction 150, as the distance 160 between the obstruction 150 and the 
wave launcher 130 increases. The graph 700 describes the probability that the inspection 
system 100 detects the obstruction 150 in a straight pipeline 140 or a curved pipeline 140. 
For example, the graph 700 represents the probability that the inspection system 100 
detects the obstruction 150 in a straight pipeline 140 when the input waveforms 235 
30 propagate at frequencies corresponding to the dominant mode. The graph 600 also 
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represents the probability that the inspection system 100 detects the obstruction 150 in a 
curved pipeline 140 when the input waveforms 235 propagate at frequencies 
corresponding to more than one mode of the input waveform 235. 

Figure 8 is a graph illustrating the probability that the inspection system 100 
5 detects a single obstruction 150 as the size (e.g., small, medium, large) of the obstruction 
150 varies. The amplitude of the reflected component 245 increases as the size of the 
obstruction 150 increases. Therefore, the probability of detection generally increases as 
the size of the obstruction 150 increases. This increase is represented by translating the 
left curve shown in Figure 8 to the right as the size of the target 150 increases. 

10 The inspection system 100 may be embodied in other specific forms without 

departing from the spirit or essential characteristics of the claimed invention. The 
foregoing embodiments are therefore to be considered in all respects illustrative rather 
than limiting on the present invention. 

What is claimed is: 
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1 1 . A pipeline inspection system comprising, 

2 a wave launcher in communication with a pipeline and adapted to transmit an 

3 input waveform having a selected input energy along a longitudinal axis of said pipeline, 

4 and to receive a reflected component of said input waveform from said pipeline, said 

5 reflected component having a characteristic reflected energy, 

6 an analyzer in communication with said wave launcher and adapted to generate 

7 said input waveform, and to receive said reflected component of said input waveform 

8 from said wave launcher, and 

9 a processor in communication with said analyzer and adapted to compare said 

10 input waveform with said reflected component of said input waveform to determine a 

1 1 characteristic of said pipeline. 

1 2. The apparatus of claim 1, wherein said processor is further adapted to compare 

2 said input waveform with said reflected component to detect an anomaly in said pipeline. 

1 3. The apparatus of claim 2, wherein said anomaly is at least one of a crack, a 

2 corrosion, a leak, a location of an end wall, an obstruction, a flange, a weld, and a 

3 restriction in said pipeline. 

1 4. The apparatus of claim 2, wherein said processor is further adapted to compare 

2 said input waveform with said reflected component to determine a location of said 

3 anomaly in said pipeline. 

1 5. The apparatus of claim 2, wherein said processor is further adapted to compare 

2 said input waveform with said reflected component to determine a shape of said anomaly 

3 in said pipeline. 

1 6. The apparatus of claim 2, wherein said processor is further adapted to compare 

2 said input waveform with said reflected component to determine one of an absolute size 

3 of said anomaly and a relative size of said anomaly relative to an internal diameter of said 

4 pipeline. 

1 7. The apparatus of claim 1, wherein said processor is further adapted to compare 

2 said input waveform with said reflected component to determine an axial curvature in 

3 said pipeline. 
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1 8. The apparatus of claim 1, wherein said processor is further adapted to compare 

2 said input waveform with said reflected component to determine location points along 

3 said pipeline relative to an initial known location. 

1 9. The apparatus of claim 1 , wherein said wave launcher further comprises a probe 

2 antenna, said probe antenna adapted for transmitting said input waveform into said 

3 pipeline. 

1 10. The apparatus of claim 9, wherein said probe antenna of said wave launcher is in 

2 physical contact with said pipeline. 

1 11. The apparatus of claim 1, wherein said analyzer is further adapted to detect said 

2 reflected component along said longitudinal axis of said pipeline. 

1 12. The apparatus of claim 1, wherein said processor is further adapted to generate a 

2 mathematical model representative of said pipeline. 

1 13. The apparatus of claim 12, wherein said mathematical model is ideal. 

1 14. The apparatus of claim 12, wherein said mathematical model is lossy. 

1 15. The apparatus of claim 12, wherein said mathematical model is one of an 

2 averaging model and a cross-sectional model. 

1 16. The apparatus of claim 12, wherein said processor is further adapted to generate a 

2 model transfer function relating a model input waveform to a model reflected component, 

3 an actual transfer function relating an actual input waveform to an actual reflected 

4 component, and to determine said characteristic at least in part by comparing said model 

5 transfer function with said actual transfer function. 

1 17. The apparatus of claim 12, wherein said processor is further adapted to determine 

2 said characteristic of said pipeline at least in part by comparing an actual reflected 

3 component with a model reflected component. 

1 1 8 The apparatus of claim 1 , wherein said analyzer is further adapted to extract a 

2 characteristic energy and phase for said input waveform and said reflected component, 

1 19. The apparatus of claim 1 , wherein said analyzer is further adapted to generate said 

2 input waveform with a frequency above a characteristic cutoff frequency of said pipeline. 

1 20. The apparatus of claim 1 , wherein said analyzer is further adapted to generate said 

2 input waveform at a frequency within a range of frequencies for which a dominant mode 

3 for said pipeline exists. 
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1 21. The apparatus of claim 20, wherein said input waveform comprises a plurality of 

2 input signals within said range of frequencies. 

1 22. The apparatus of claim 21 , wherein said analyzer is further adapted to detect 

2 differences in velocity between said plurality of input signals as said input signals 

3 propagate in said pipeline, and said processor is further adapted to determine a curvature 

4 of said pipe along said longitudinal axis from said differences in velocity. 

1 23. The apparatus of claim 21 , wherein said analyzer is further adapted to detect 

2 differences in velocity between reflected components of each of said plurality of input 

3 signals to determine a curvature of said pipeline along said longitudinal axis. 

1 24. The apparatus of claim 1 , wherein said analyzer is further adapted to generate an 

2 electromagnetic waveform as said input waveform. 

1 25. The apparatus of claim 1 , wherein said analyzer is further adapted to generate an 

2 acoustic waveform as said input waveform. 

1 26. The apparatus of claim 1 , wherein said analyzer is further adapted to generate said 

2 input waveform as one of a spread spectrum waveform, a chirp waveform, and a soliton 

3 waveform. 

1 27. The apparatus of claim 1, wherein said analyzer is further adapted to generate said 

2 input waveform as a wideband waveform. 

1 28. The apparatus of claim 1 further comprising calibration elements adapted to 

2 temperature stabilize said analyzer. 

1 29. A method of inspecting a characteristic of a pipeline, said method comprising, 

2 transmitting an input waveform having a selected input energy along a 

3 longitudinal axis of said pipeline, 

4 receiving a reflected component of said input waveform from said pipeline, said 

5 reflected component having a characteristic reflected energy, and 

6 comparing said input waveform with said reflected component of said input 

7 waveform to determine said characteristic of said pipeline. 

1 30. The method of claim 29 further comprising, comparing said input waveform with 

2 said reflected component to detect an anomaly in said pipeline. 
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1 31. The method of claim 30, wherein said anomaly is at least one of a crack, a 

2 corrosion, a leak, a location of an end wall, an obstruction, a flange, a weld, and a 

3 restriction in said pipeline. 

1 32. The method of claim 30 further comprising, comparing said input waveform with 

2 said reflected component to determine a location of said anomaly in said pipeline. 

1 33 . The method of claim 30 further comprising, comparing said input waveform with 

2 said reflected component to determine a shape of said anomaly in said pipeline. 

1 34. The method of claim 30 further comprising, comparing said input waveform with 

2 said reflected component to determine one of an absolute size of said anomaly and a 

3 relative size of said anomaly relative to an internal diameter of said pipeline. 

1 35. The method of claim 29 further comprising, comparing said input waveform with 

2 said reflected component to determine an axial curvature in said pipeline. 

1 36. The method of claim 29 further comprising, comparing said input waveform with 

2 said reflected component to determine location points along said pipeline relative to an 

3 initial known location. 

1 37. The method of claim 29, further comprising, detecting said reflected component 

2 along said longitudinal axis of said pipeline. 

1 38. The method of claim 29 further comprising, generating a mathematical model 

2 representative of said pipeline. 

l 39. The method of claim 38, wherein said mathematical model is ideal. 

1 40. The method of claim 38, wherein said mathematical model is lossy. 

1 41 . The method of claim 38, wherein said mathematical model is one of an averaging 

2 model and a cross-sectional model. 

1 42. The method of claim 38 further comprising, generating a model transfer function 

2 relating a model input waveform to a model reflected component, an actual transfer 

3 function relating an actual input waveform to an actual reflected component, and to 

4 determine said characteristic at least in part by comparing said model transfer function 

5 with said actual transfer function. 

1 43. The method of claim 38 further comprising, determining said characteristic of 

2 said pipeline at least in part by comparing an actual reflected component with a model 

3 reflected component. 
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1 44 The method of claim 29 further comprising, extracting a characteristic energy and 

2 phase for said input waveform and said reflected component. 

1 45. The method of claim 29 further comprising, generating said input waveform with 

2 a frequency above a characteristic cutoff frequency of said pipeline. 

1 46. The method of claim 29 further comprising, generating said input waveform at a 

2 frequency within a range of frequencies for which a dominant mode for said pipeline 

3 exists. 

1 47. The method of claim 46, wherein said input waveform comprises a plurality of 

2 input signals within said range of frequencies. 

1 48. The method of claim 47 further comprising, detecting differences in velocity 

2 between said plurality of input signals as said input signals propagate in said pipeline, 

3 and determining a curvature of said pipe along said longitudinal axis from said 

4 differences in velocity. 

1 49. The method of claim 47 further comprising, detecting differences in velocity 

2 between reflected components of each of said plurality of input signals to determine a 

3 curvature of said pipeline along said longitudinal axis. 

1 50. The method of 29 further comprising, generating an electromagnetic waveform as 

2 said input waveform. 

1 51. The method of claim 29 further comprising, generating an acoustic waveform as 

2 said input waveform. 

1 52. The method of claim 29 further comprising, generating said input waveform as 

2 one of a spread spectrum waveform, a chirp waveform, and a soliton waveform. 

1 53. The method of claim 29 further comprising, generating said input waveform as a 

2 wideband waveform. 

1 54. The method of claim 29 further comprising, calibrating said analyzer to be 

2 temperature stable. 

1 55. A method of determining a location of a point along a pipeline, said method 

2 comprising, 

3 transmitting an input waveform having a selected input energy along a 

4 longitudinal axis of said pipeline, 
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5 receiving a reflected component of said input waveform from said pipeline, said 

6 reflected component having a characteristic reflected energy, and 

7 comparing said input waveform with said reflected component of said input 

8 waveform to determine said location of said point along said pipeline. 

1 56. A method of inspecting a characteristic of a pipeline, said method comprising, 

2 transmitting an input waveform having a selected input energy along a 

3 longitudinal axis of said pipeline, 

4 receiving a reflected component of said input waveform from said pipeline, said 

5 reflected component having a characteristic reflected energy, and 

6 determining said characteristic of said pipeline using an error estimate, said error 

7 estimate depending on a known point along said pipeline relative to said characteristic. 

1 57. The method of claim 56 wherein said error estimate further depends on one of a 

2 conductance, a radius, and a cross-sectional shape of said pipeline. 

1 58. A method of inspecting a characteristic of a pipeline, said method comprising, 

2 generating an input waveform, 

3 launching said input waveform into said pipeline, 

4 receiving from said pipeline a reflected component having a characteristic 

5 reflected energy of said input waveform, 

6 calculating a mathematical function of said characteristic reflected energy from 

7 said reflected component of said input waveform, 

8 determining a model mathematical function of model reflected energy from a 

9 model component of a model input waveform, and 

10 determining said characteristic of said pipeline by comparing said mathematical 

1 1 function of said reflected energy to said model mathematical function of said model 

12 reflected energy. 
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Abstract of the Invention 

The invention is directed to a system and method for non-invasive pipeline 
inspection. According to one embodiment, the system includes a processor, an analyzer, 
and a wave launcher. The wave launcher is adapted to transmit an input wideband 

5 waveform having a selected input energy into the pipeline along a longitudinal axis, and 
to receive from the pipeline a reflected component of the input waveform having a 
reflected energy. The analyzer is adapted to generate the input waveform, and to receive 
the reflected component of the input waveform from the wave launcher. The processor is 
adapted to compare the input waveform with the reflected component of the input 

10 waveform to determine characteristics. 
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P.O. Address (line 2) 


City 


Ashland 


State 


MA 


ZIP 


01721 


Country 


USA 


Q Additional inventors are being named on the 


supplemental Additional Inventor(s) sheet(s) attached hereto. 


Name of Additional Joint Inventor, if any: 


□ A petition has been filed for this unsigned inventor 


Given Name (first and middle [if any]) 


Family Name or Surname 


Claude P. 


Brancart 


Inventor's Signature 




Date 




Residence 


City 


Brunswick 


State 


ME 


Country 


USA 


Citizenship 


US 


Post Office Address 


18 Juniper Road 


P.O. Address (line 2) 


City 


Brunswick 


State 


ME 


ZIP 


04011 


Country 


USA 


Name of Additional Joint Inventor, if any: 


D A petition has been filed for this unsigned inventor 


Given Name (first and middle [if any]) 


Family Name or Surname 






Inventor's Signature 




Date 




Residence 


City 




State 




Country 




Citizenship 




Post Office Address 




P.O. Address (line 2) 


City 




State 




ZIP 




Country 





BIANCOJV6209/7.1046587-1 



